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MEASUREMENT  OF  PROPERTIES  OF  SPREAD  CHANNELS 
BY  THE  TWO- FREQUENCY  METHOD 
WITH  APPLICATION  TO  RADAR  ASTRONOMY 


I.  INTRODUCTION 

In  many  communication  channels,  even  in  the  absence  of  additive  noise,  a  signal  applied  to 
the  channel  might  be  seriously  distorted.  This  distortion  arises  because  the  channel  destroys 
the  original  phase  and  amplitude  relationships  between  the  various  components  of  the  input  signal. 
If  a  detailed  knowledge  of  the  phase  and  amplitude  characteristics  of  the  channel  exists,  the  dis¬ 
tortion  may  be  removed  by  appropriate  filtering.  In  addition  to  the  distortion,  there  will  also  be 
a  certain  amount  of  time  variation  in  many  channels.  The  application  of  a  single  sine  wave  to 
the  channel  will  therefore  result  in  an  output  signal  of  finite  bandwidth.  Because  of  these  prop¬ 
erties,  it  is  most  convenient  to  describe  such  a  channel  by  certain  statistical  properties.  In 
this  report  we  shall  assume  throughout  that  these  statistical  properties  are  stationary  in  time. 

As  a  result,  the  various  different  frequency  components  in  the  output  will  be  completely  uncor¬ 
related  when  the  input  waveform  is  a  single  sine  wave.  Similarly,  it  will  be  assumed  in  the 
analysis,  wherever  convenient,  that  the  statistical  properties  of  the  channel  are  largely  inde¬ 
pendent  of  the  radio  frequency  of  observation.  This  is  equivalent  to  assuming  that  the  signals 
scattered  at  various  delays  are  uncorrelated,  and  this  is  often  a  very  good  approximation. 
However,  such  an  assumption  is  not  necessary  for  the  application  of  the  measuring  scheme. 

There  is  considerable  literature  on  the  description  and  measurement  of  such  doubly  spread 
1 

channels.  Green  describes  such  a  channel  in  terms  of  a  target-scattering  function  which  de¬ 
scribes  the  way  in  which  the  power  incident  is  distributed  in  delay  and  doppler  displacement. 
Kailath  characterized  this  type  channel  in  terms  of  a  tap-gain  correlation  function  which, 

loosely  stated,  describes  the  autocorrelation  function  of  the  time-varying  impulse  response  as 

3 

a  function  of  delay.  Hagfors  describes  the  channel  in  terms  of  a  two-dimensional  correlation 
function  describing  the  loss  of  correlation  between  two  time-varying  frequency  responses  with 
a  shift  in  frequency  Af  and  a  shift  in  time  T.  An  instructive  survey  of  the  various  descriptions 

4 

has  recently  been  given  by  Gallager. 

Several  methods  have  been  proposed  for  the  measurement  of  the  various  functions  describing 

the  statistical  properties  of  the  spread  channel.  Kailath  proposes  to  use  a  noise  waveform  at 

the  input  and  compute  the  fourth -order  moments  at  the  output  to  derive  the  tap-gain  correlation 

function.  This  method  appears  to  work  both  when  the  product  of  doppler  spread  B  and  delay 

1  3 

spread  L  is  greater  than  and  less  than  unity.  Green  and  Hagfors  both  describe  methods  of 
measurement  that  only  work  well  for  underspread  channels,  i.e.,  when  BL  <  1;  however,  they 
did  not  consider  the  presence  of  additive  noise.  Price^  and  Levin7  dealt  with  the  problem  of 
additive  noise  in  the  estimation  of  the  statistical  properties  of  the  channel. 


1 


In  this  report,  we  analyze  a  particular  measuring  scheme  based  on  a  two-frequency  trans¬ 
mission  when  additive  noise  is  present.  It  will  be  shown  that  the  two-frequency  method  of  ob- 

5 

servation,  apparently  contrary  to  common  belief,  can  be  made  to  work  satisfactorily  both  in  the 

overspread  and  underspread  situations.  A  similar  discussion  has  recently  been  carried  out  for 

4 

the  case  of  pseudo-random  or  a  chirp  input  to  the  random  channel. 

H.  DESCRIPTION  OF  MEASUREMENT  SCHEME 

The  basic  observation  scheme  is  shown  in  block  diagram  form  in  Fig.  1.  The  transmitter 
emits  two  sine  waves  of  equal  strength  at  frequencies  fQ  ±  Af/2.  The  modulation  frequency  Af 
must  be  variable  and  it  is  highly  desirable  that  Af  is  very  accurately  known,  particularly  when 
there  are  large  delays  in  the  channel  such  as  in  planetary  and  lunar  radar  astronomy.  It  is 
convenient  to  choose  Af  as  an  integral  multiple  of  some  smallest  frequency  separation  AfQ,  and 
let  this  be  derived  from  a  frequency  standard.  The  necessity  of  this  requirement  will  be  briefly 
explained  in  more  detail  below. 

The  received  signal  will  consist  of  two  spectra  which,  in  the  absence  of  mean  doppler  dis¬ 
placement  which  we  will  ignore  in  this  report,  will  be  centered  on  the  two  transmitted  frequencies 
fQ  ±  Af/2.  If  a  mean  doppler  displacement  is  present  we  shall  assume  that  the  local  oscillators 


Fig.  1.  Block  diagram  of  two-frequency  measurement  scheme. 
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of  the  receiving  system  are  programmed  to  correct  for  it.  When  the  channel  is  "underspread," 
nearly  all  the  significant  measurements  can  be  made  without  the  two  spectra  overlapping;  in  the 
"overspread”  situation  this  will  not  be  the  case.1  The  received  signal  is  multiplied  by  a  distorted 
replica  of  the  transmitted  signal,  the  distortion  amounting  to  a  change  in  the  modulation  frequency 
by  an  amount  f^,  so  that  the  modulation  frequency  of  the  multiplying  signal  is  fL  +  Af.  Also,  the 
center  frequency  of  the  multiplying  signal  is  different  from  fQ  by  an  amount  f^  =  fQ  —  f  1  •  The 
offset  frequency  f ^  must  be  chosen  so  that  it  is  greater  than  the  maximum  bandwidth  Bm  of  the 
response  of  the  time-varying  channel  when  a  single  sine  wave  is  applied  at  the  input.  The  IF 
frequency  f^  must  be  greater  than  the  sum  of  f ]V2,  the  maximum  frequency  separation  Af,  and 
half  the  maximum  bandwidth  in  order  to  prevent  "foldover"  (see  Fig.  2).  Following  the  multiplier 
the  signal  is  split  into  two  parts  by  upper-  and  lower-band  filters.  In  the  strictly  under  spread 
situation,  these  two  filters  should  be  constructed  so  that  only  the  signals  centered  on  f^  +  fjV2 
and  f^  —  fjV2  are  passed,  respectively.  In  the  overspread  situation  such  a  discrimination  will 
no  longer  be  possible,  and  we  must  allow  the  filters  also  to  pass  some  of  the  signals  centered 
on  f^  ±  [(L/2)  +  Af].  This  will  tend  to  introduce  a  certain  amount  of  "self-noise"  into  the 
system. 

The  upper-band  signal  is  passed  into  a  delay  line  with  discrete  delays  0,  T^,  T^,  etc., 

corresponding  to  the  time  shifts  in  the  two-dimensional  correlation  function.  These  time  shifts 

.  \ 

should  be  chosen  to  be  integral  multiples  of  fj^  for  convenience.  The  delayed  and  undelayed 
signals  are  multiplied  together  and  the  amplitude  and  phase  of  the  component  in  the  output  cor¬ 
responding  to  the  frequency  fT  are  determined  by  means  of  the  operations  shown  in  Fig.  1.  The 

^  -1 
integration  time  0  should  be  chosen  as  an  integral  multiple  of  the  period  AfQ  of  the  lowest 

modulation  frequency  for  reasons  that  will  become  apparent  from  the  analysis  that  follows. 

The  pairs  of  numbers,  or  the  set  of  complex  numbers,  resulting  from  the  measurement  are 
shown  below  to  result  in  a  determination  of  the  complex  correlation  functions  with  frequency 
shift  Af  and  with  time  shifts  0,  T^,  T^,  T^,  etc. 

Let  us  now  briefly  return  to  the  reason  why  the  modulation  frequency  must  be  so  accurately 
known.  Consider  as  an  example  the  situation  encountered  in  radar  astronomy.  The  transmission 
of  the  two  sine  waves  corresponds  to  measuring  one  particular  Fourier  component  of  the  target¬ 
scattering  function,  giving  power  as  a  function  of  range.  If  this  function  is  to  be  recovered 
completely  we  must  know  both  amplitude  and  phase  of  the  Fourier  components  for  different  Af. 


Fig.  2.  Various  frequency  spectra 
in  the  processing  scheme. 


1 3-31-30141 


(a) 


(b) 


(c) 


TRANSMITTED  SPECTRUM 


'0 

-  Af  - 


A  A 

fo  ^  kF 


RECEIVED  SPECTRUM 


h - ^  ■ 


-  Af  - H 


A  A  .  A  A 

-H  Bm  k-  f,  FREQUENCY 


SPECTRUM  AFTER  MULTIPLICATION 
BY  DISTORTED  REPLICA 


3 


The  phase  and  amplitude  relationships  must  be  well  defined  at  the  target.  In  radar  astronomy 

and  in  other  applications  there  will,  however,  be  a  large  number  of  modulation  cycles  between 

the  transmitter  and  the  target.  A  very  slight  inaccuracy  in  the  modulation  frequency  might, 

therefore,  result  in  a  large  phase  inaccuracy  at  the  target. 

The  sets  of  results  obtained  with  a  number  of  different  frequency  separations  Af  can,  if 

3 

desired,  be  used  to  compute  the  target-scattering  function  by  a  double  Fourier  transform 
(see  Sec.  III-E  below). 


III.  ANALYSIS  OF  MEASURING  SCHEME 


In  this  section  it  will  first  be  shown  how  the  mean  values  of  the  output  of  the  processing 
scheme  shown  in  Fig.  1  are  related  to  the  two-frequency  correlation  function  both  for  under¬ 
spread  and  overspread  channels,  first  without  additive  thermal  noise,  then  with  thermal  noise 
included.  In  the  subsequent  sections  we  will  repeat  the  analysis  under  the  same  situations  for 
the  variance  of  the  output  so  that  confidence  limits  can  be  placed  on  the  mean  values  determined. 

Before  starting  on  this  program  it  is  worthwhile  to  define  the  complex  two-frequency  cor¬ 
relation  function  in  terms  of  a  time-varying  frequency  response.  Suppose  the  complex  amplitude 
of  a  pure  sine  wave  at  the  channel  input  is  e ^  and  that  the  frequency  is  fQ.  The  complex  amplitude 
at  the  output  is  then  time  varying  because  of  the  channel  fluctuation  and  can  be  expressed  as 
follows: 


e  ,(t)  =  H(fn,  t)  e. 
out  0'  in 


(1) 


The  ratio  of  mean  power  at  the  output  to  input  power  becomes 

IT*  =  <|H(f0,t)|2>  .  (2) 

in  & 

where  the  averaging  is  considered  an  ensemble  averaging,  in  general.  The  (unnormalized)  two- 
frequency  correlation  function  is  now  defined  as 

R(Af,  T)  =  <H*(f0  -  ^,t)  H(f0  +  +  T)>avg  .  (3) 

where  we  have  assumed  that  the  channel  is  stationary  both  in  time  and  frequency.  In  case  the 
channel  is  not  stationary  in  frequency,  the  correlation  function  will  also  depend  on  fQ.  The 
transmission  loss  is  obtained  from  Eq.  (3)  by  substituting  Af  =  0  and  T  =  0.  We  assume  through¬ 
out  this  report  that  H(f^,  t)  is  a  Gaussian  process. 

A.  Evaluation  of  Mean,  No  Noise 

The  signal  at  the  input  of  the  receiver  will  have  the  complex  amplitude  referred  to  the 
center  frequency  fQ  given  by 


F.  (t)  =  H(fn 
in  0 


Af  .  -i7rAft 
e 


+  H(f0 


,  Af  .  in-Aft 

+  ~r>  t)  e 


(4) 


provided  the  transmitted  amplitudes  are  equal  to  unity.  The  complex  amplitude  at  the  output  of 
the  first  multiplier,  this  time  referred  to  the  intermediate  frequency  f^  =  f^  —  f^,  becomes 
(apart  from  a  constant  of  proportionality) 


4 


F2(t)  =  [H(f0  -  t)  +  H(f0  +  t)  e2,rlAft] 

+  [H(f0  +  ^,  t)  +  H(f0  -  ^ ,  t)  e'2,riAft]  e  ^  ljt 


(5) 


The  upper-  and  lower-band  signals  appearing  separately  above  can  always  be  separated  in  the 
two  filters  following  the  first  multiplier,  leaving  two  signals  with  amplitudes  referred  to 
frequencies  f^  +  fjV2  and  *2  respectively: 


F2  (t)  =  [H(f0  -  t)  +  H(f0  +  t)  e2,riAft] 


F‘(t)  =  [H(f0  +  ^,  t)  +  H(f0  -  t)  e"2,riAftJ 


(6) 


We  next  filter  out  as  much  of  the  undesired  signal  as  possible  by  assuming  the  maximum  band¬ 
width  of  the  time-varying  response  to  be  B  ,  so  that  the  following  representation  applies. 

H(f0  ±  t)  =  J  H(f0  ±  v)  e2*1"1  dv 


fBm/2  H(fn  ±  ¥-.v)  e2nil,t 


dv 


(7) 


By  including  rectangular  filters  in  the  lower-  and  upper-band  filters  of  bandwidth  equal  to  B  , 
the  desired  signal  is  passed  undistorted  and  a  certain  fraction  of  the  undesired  signal  is  rejected. 
The  results  of  passing  the  signals  [Eq.  (6)]  through  such  filters  centered  on  f^  ±  f^/2  are: 


F3+(t)  =  H(f0  -  ^.t)  + 


rBn/2 

J-(B  /2)+Af 
m 


H(f„  +  4^.  v  -  Af>  eZnltv  dv  Af  <  B 
U  c.  ’  m 


Af  >  B 


(8a) 


Af 


F3(t)  =  H(f0  +  ^.t)  + 


(MB  /2)-Af  A,  ,  .. 

J  m  .  H(fn  -¥■.»*  A«  e2lrlt"  d,  Af  <  B_ 


-(Bm/2) 


Af  >  B 


(8b) 


The  second  multiplication  with  subsequent  coherent  detection  of  the  component  at  frequency  fj 
is  represented  mathematically  as: 


r6 

V,T<> =  X 


F3  *  (t)  F*  (t  +  Ty)  dt 


(9) 


or  explicitly. 


dt  |H*(fn  -  4^.t)  +  D(Bm  -  Af)  \  m'  H*(f_  +  ^,v  -Af)  e 

02  m  J-(B  /2)+Af  °  2 

m 


» ■  £ 

X  JhUq  +  t  +  Tf)  +  D(Bm  -  Af) 


IT */r  ,  Af  ..  A  r*\  _  —  27ril^ t 


d^ 


(Bm/2)-Af 

(B  /  2) 
nr 


+  Af)  e 


27ri^(t+Tj 


dv 


(10) 
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where  D(x)  is  a  function  that  is  zero  when  the  argument  is  negative,  and  unity  when  it  is  positive. 
The  mean  value  of  the  output  becomes 

<V(Tl)>avg  =  6  <H*(fo"  TT*1)  H<fo  +  ;T''t  +Ti)>avg  =  eR(Af.T|>  •  dD 

The  other  three  terms  will  vanish  exactly  if  the  integration  time  0  is  chosen  as  an  integral 
multiple  of  Af  ” *,  as  explained  in  Sec.  II.  This  can  be  seen  from  the  following  development. 


<r 


(B  /2)-Af 


Af  2iriy(t+T  )  \ 

H(f-  +  Af)  e  dMdt  ) 

(B_/2)  °  2  / 


avg 


\  dtLdv,y 


r(B  /2)-Af 


m 


Af 


Af 


<H*(f0  -■=£,*')  H(fn  —  v  +  Af)> 


-27ri[y 1  — ^ )  t— T#  ] 


0  2 


avg 


dy 


r(B  /2)-Af  Af  7  ZitivT. 


avg 


-2?riAft 
e  dt 


=  0 


(12) 


and  the  latter  of  these  two  integrals  obviously  must  vanish  exactly  under  the  conditions  imposed 
on  0.  The  product  of  the  two  integrals  in  Eq.  (10)  will  also  vanish  for  the  same  reason.  There 
is  only  one  exception  to  this,  namely,  when  Af  =  0  all  four  terms  in  Eq.  (10)  contribute  an  equal 
amount.  Hence,  we  conclude. 


<V(Af.Ti)>avg=<v(mAf0, 


8,K'() 

m  ^  1  ,  l  =  0,  1,  .  .  . 

4eR(o.ji) 

cT 

ii 

o 

It 

£ 

(13) 

These  results  apply  whether  or  not  the  channel  is  overspread. 


B.  Evaluation  of  Mean  with  Noise  Present 

When  additive  noise  is  present  at  the  receiver  input  there  will  be  an  additional  term  in  the 
input  signal. 

F.n(t)  =  H(f0  -  ^ ,  t)  e-,riAft  +  H(f0  +  ^ ,  t)  effiAft  +  n(t)  (14) 

Here,  n(t)  is  a  noise  signal  which  might  be  expressed  as 


rfn+(Af  /2)  ,  ... 

n(t)  =  \  maX  N<f>  e  df  -  H5) 

Jfn-(Af  /  2) 

0  max7 

where  Af max  signifies  a  largest  modulation  frequency  to  be  used.  In  the  output  of  the  processor 
certain  new  noise  terms  will  appear.  The  terms  containing  cross  products  of  noise  and  signal 
will  not  contribute  to  the  mean,  because  noise  and  signal  are  uncorrelated.  The  noise-noise 
term,  however,  might  contribute  to  the  mean  output.  To  see  this  we  may  argue  as  follows. 

The  noise  associated  with  F^(t)  becomes 
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psX  - 


n  -(Af/2)+(Bm/2) 
^-(Af/2)-(Bm/2) 


N(v  +  f^)  e2lril,t  dt> 


(16a) 


and,  similarly,  for  F^(t), 


r(Af/2)+(B  / 2)  ,  . 

F3_(t)n  =  V  m  N(t>  +  f  )  e2,rll't  du  . 

J(Af/2)-(Bm/2)  0 

The  noise  contribution  to  the  outputsignal  V(T^),  i.e.,  Vn(Tj),  becomes 

■  r d*  i 


9  r-(Af/2)t(Bm/2| 

N*(y  +  fn)  e 

(Af/2)-(Bm/2) 


C 

Jt 


(Af/2)+(B  /2) 


'(Af/2)-(Bm/2) 


+  fj  e 


27rii',(t+Ti) 


(16b) 


(17) 


Taking  the  mean  of  this  we  find 


<Vn<T,»a»g-£  *£ 


0  r-(Af/2)+(B  / 2)  ,  ,  .  . 

r  m  d„  <|N(,  +f0)|2>  e’2™1 


(Af/2)-(Bm/2) 


r(Af/2)+(B  /2)  2irit',(t+T.) 

X  \  dr'  e  6(i>'  -  v) 

J(Af/2)-(B  /2) 


0  r-(Af/2)+(Bm/2) 

(Af/2)-(B  /2) 


<|  N(i>  +  f  )  |  >  <\v  e 


27ri^T. 


X  {D  [If  -  J  (Af  -  BJ]  -  D  [v  -  j  (Af  +  B  ) ]} 


sin7r  T#(B  —  Af) 

(B  —  Af)  0N  D(B  —  Af)  -m7Tf(— 

m  o  m  7t1  .  ( B —  Ai) 

£  m 


where  Nq  is  the  noise  power  per  unit  bandwidth,  and  D(x)  has  the  same  significance  as  in  Eq.  (10). 
Hence,  the  output  of  the  processor  with  noise  present  becomes 


<V(mi,0'  Xav,,  *  9 


RHf0*  ff)  +  No^m-^V  - X 


sin7r  fT  (Bm"mAf0) 

Li 


,  (B  -mAL) 
fL  m  0 


D(Bm-mAf0) 


m  ^1  (18) 

We  conclude  that  the  processor  gives  an  output  proportional  to  the  desired  correlation  func¬ 
tion  whenever  the  two  lines  can  be  completely  separated.  When  the  two  lines  overlap  there  will 
be  a  noise  bias  in  addition  (see  Fig.  3);  hence,  this  bias  must  be  established  before  the  correla¬ 
tion  function  can  be  determined. 
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1 3-3 1->0I5 1 


Fig.  3.  Output  of  the  processor  shown 
in  Fig.  1  when  noise  is  present. 


C.  Evaluation  of  Variance,  No  Noise 

The  variance  of  the  output  may  be  computed  directly  from  the  expression  [Eq.  (9)]: 

<lV(Ti)|2>avg  =  ^  ^  dtdt'  <F3*(t)  F3"(t  +  Ti>  F3+(t,)  F3'"(t’  +  Ti»  avg  '  <19> 

Now,  since  we  are  assuming  a  Gaussian  process,  the  integrand  can  be  expanded  in  terms  of 
second-order  moments  as  follows. 

<F3+*(t)  F3'(t  +  Tt)  F3(t')  F^’*(tl  +  Tf)>avg 

=  <F^*(t)  F3-(t  +  Ti»avg  <F3+(t')  F3-*(f  +  avg 
+  <F3+*(t)  F3+(F)>avg  <F;(t  +  T,)  F;*(f  +  Ti»avg 

+  <F^(t)  F3'*(f  +  T,)> avg  <F3"(t  +  T,)  F3+<t')>avg  .  (20) 

The  first  term  on  the  right-hand  side  in  Eq.  (20)  can  be  seen  to  correspond  simply  to  the  square 
of  the  mean  value  which  has  already  been  determined.  The  last  term  in  Eq.  (20)  is  completely 
negligible  because  it  involves  mean  values  of  products  of  frequency  responses,  either  both 
unconjugated  or  both  conjugated.  Only  the  middle  term  will  therefore  contribute  to  the  variance. 
In  order  to  evaluate  this  term  we  write 

g  /2 

F+(t)  =  H(f  -  ^r,t)  +  D(B  -  Af)  f  m  H(f.  +  4^,  v  -Af)  eZvivt  dv  , 

3  u  2  m  J_(B  /2)+Af  o  2 

m 

F.T(t)  =  H(fn  +  Q?-,  t)  +  D(B  —  Af)  C<Bm//  f  H(fn  -  i/  +  Af)  e2™*  dr  ,  (21) 

3  02  m  J-(B  /  2)  0  2 

m 

where  the  factor  D(B  —  Af)  is  unity  whenever  B  >  Af,  and  zero  otherwise.  We  have: 
m  m 
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<F**(t)  F+(t')>  „  =  R(0,t'  -t)  +  D(B  -Af)  CB-/2  di/'  C  111  dP 

3  3  aVg  m  J-<Bm/2)  J-(Bm/2)+Af 


-B  /2 
nr 


X  e'2’rl(‘''t‘l't')  <H*(f0  H(fn  +  ^-.v  —  Af)> 


+  D(B  —  Af) 
m 


I 


B  /2 


m 

l  , 
m 


dy*  \  *“  dv  e 

(B/Z)  J-(B  /2)+Af 


B  /2 
nv 


avg 
27ri(i'  't'  —  t>t) 


B  /2 


x<H(f0-^[-‘',>H*<fo  +  T '  "  -  Af>avg  +  D<Bm  -  Af)  H  d"dl 


-(Bm/2)+Af 


x<H^f0  +  f ,  ,-Af)  H(f0  +  f .  v '  —  Af  )>  avg  e 
Because  of  the  assumed  time  stationarity  of  the  statistical  properties  we  must  have 


Af 


iZ7ri(i/,t,-^t) 


(22) 


<H*<f0  -  Tf  -  H<f0  +  "^avg  =  <H*<f0  -  ¥•  v)  H<f0  +  ¥'  "»avg  6(*"  " 


Using  this  we  obtain 


<F3+*(t)  F3+{t’)>avg  =  R(0,F  -t)  +  D(Bm-Af) 


£ 


B  /2 

rrv 


(B  /2)+Af 


du  e 


27ri^(t'  -t) 


X  <H*(fQ  -  v  -  Af)  H(f0  +  v  -  Af>av  e2,riAft  +  D(Bm  -  Af)  e_27riAft' 


£ 


B  /2 
nr 

<Bm/2)+Af 


e2»iv(t'-t)  <H*(fo  +  _Af)  H(f0  -  -Af)>  dp 


+  D(B  -  Af) 


rB  /2 
(  m 


Av  gZjTH'tt'-t)  <H*(f  „  —  Af)  H(f„,  v  —  Af)>  .  (22a) 


m  '  J  (B  /2)+Af  '  °'  °*  ayg 


m 


The  other  moment  which  must  be  evaluated  is 


<F"(t  +  Tf)  F”  (t1  +  Ti)>  =  <F3"(t)  F"  (t’)> 


avg 


(23) 


The  last  equation  follows  from  the  time  stationarity  of  the  random  process.  The  average  in 
Eq.  (23)  can  be  obtained  from  Eq.  (22)  by  changing  the  sign  on  Af  in  all  arguments  of  H,  by 
interchanging  the  roles  of  t  and  t1,  and  by  changing  the  limits  of  integration.  One  then  obtains 


<F"(t)  F"*(t')>avg  =  R(0,  t  —  t 1 )  +  D(Bm-Af)  e'2,rlAft 


■£ 


(Bm/2)-Af  d(/  e-2»U»(f-t) 

(B  /2) 


X  <H*(f0  +  v  +  Af)  H(f0  -  v  +  Af)>ayg  +  D(Bm  -  Af)  e2,rlAft  ^  cb 


Af 


-(Bj2) 


X  e 


-  27rii^  ( t1  -t) 


<H*(f0  -%.V+  Af)  H(f0  +  +  Af)>avg  +  D(Bm  -  Af)  J 


(Bm/2)-Af 

(B  /  2) 
m 


dv 


x  e 


—  27r  ii^  (t 1  -t) 


<1  H(fQ,  v  +  Af)  |  2> 


avg 


(24) 
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Of  the  16  terms  resulting  from  the  multiplication  of  Eq.  (22)  by  Eq.  (24),  only  those  terms  which 
depend  on  the  combination  t  —  t*  will  contribute  in  the  integrals  over  t  and  t'.  The  terms  depend¬ 
ing  explicitly  on  either  t  or  t1  will  all  vanish.  The  contributing  terms  will  be  the  following  six: 


(I)  |R(0,t  —  t')  |  2  + 

p(B  / 2) -Af  „  ....  , 

(n)  +R(0,t'  —  t)  •  D(B  -Af)  l  m  du  e'2,rl't  -t)l'  <  j  H(fn,  v  +  Af)  |  2>  + 

m  J-(B  /2)  0  avg 

nr 

(HI)  +  R(0,  t  —  t1)  D(B  —  Af)  f  d!/  e2lrl(t'"t)l'  <|H(fn,  >/  -  Af)|2>  + 

J-{Bm/2)+Af  °  avg 

(IV)  +  D(B  -Af)  fBm/2  d„ 

J-(B  /2)+Af 
m 

X  <H*(f0  -  f  -  Af)  H(f0  +  -  Af»avg 

X  1(^/2)  Af  ^  e"2’ril/(t,'t>  <H^f0  +  ¥■  V  +  Af)  H(f0  "  ¥•  V  +  Af»  avg  + 


(V)  +  D(B  -Af)  fBm/2  d„  ^Mf-t) 

m  J-(Bm/2)+Af 

X  <H*(f0  +  -  Af)  H(f0  -  v  -  Af»avg 

X  l  ‘  m/  ’  dt-  e_2,rll'(t,‘t)  <H*(f  -  ^ ,  v  +  Af)  H(f  +  v  +  Af)>  + 

J.(B  /2)  0  2  0  2  avg 

m' 

r(B  /2)-Af  .  ....  ..  , 

(VI)  +D(B  -Af)\  m  di-  e'Zn,t <|H(f  v+Af)|2> 

J-(B  /  2)  0  avg 

m 

X  fBm/2  e2H(t'-t)r  <|H(f  „-Af)|2>  d^  . 

J-(B  /2)+Af  U  a  g 


(25) 


The  interpretation  and  discussion  of  such  a  complicated  expression  is  rather  difficult,  and 
it  simplifies  matters  greatly  if  we  can  make  a  specific  assumption  about  the  channel  properties 
to  insert  into  the  general  terms  in  (25).  As  a  specific  example  we  choose 


<H*(f0-^(>,)  H(f0  + 


avg 


j  TT  2  “if  Af  ifM<- 
° 

1 0  otherwise 


(26) 


This  corresponds  to  assuming  an  exponential  power-vs-range  variation  with  power  falling  off  to 
e’1  at  t  =  tq.  For  R(Af,  T),  we  obtain 


R(Af,  T)  -  4  +  ZniT  Af  J 


r  B  /2 

(  m 


(B  /2) 
m 


ZnivT  . 
e  d^  = 


oB 


sin7rTB 


m 


1  +  27riT  Af 
o 


ttTB_ 


(27) 


10 


The  contribution  to  the  variance  caused  by  the  terms  in  (25)  in  this  case  is  found  by  substitution 

of  Eqs.  (26)  or  (27)  and  integration  over  t  and  t'  between  the  limits  0  and  ©.  On  the  assumption 

that  ©  »  B  ”  ,  the  contributions  become 
m 


Term  (I): 

Terms  (II)  &  (in): 


©B  a' 
m 


2©(B  -  Af)  olZ  if  | Af  |  <  B 


m 


otherwise 


Terms  (IV)  &  (V): 


2©(B  —  Af)  a1 


m 


1  +  (27rTQAf) 


2  if  |  Af  |  <  B 


m 


otherwise 


Term  (VI): 


e(Bm-Af)a2  if|Af|<B 


otherwise  .  (28) 

Half  this  variance  corresponds  to  the  real  and  half  to  the  imaginary  part  of  R(Af,  T),  so  we 
might  define  the  relative  uncertainty  of  the  two  components  as  follows  (Af  >  0). 

,1/2 


<5 (Re  V) 


<Re V> 


avg 


^1/2  jl  +  (2irroAf)2  +  2D(Bm  -  Af)  (l  -  [3  +  { 2tt tqA f) 2 ] 

/sin7rTB  \ 

7Bme(  ?rTBm  ) 


g(ImV) 

(imV) 

avg 


^1/2  Jl  +  (27TToAf)2  +  2D(Bm  -  Af)  (l  -  [3  +  U^Af)2]  J^ 


1/2 


/sin7rTB 

I  r 


ttTB 


(29) 


(27TToAf) 


For  the  singular  case  of  Af  =  0,  one  obtains 


6(Re  V) 


<Re  V> 


avg 


, _  /Sin 7rTB  v 

7B mG  (  ttTB  ) 

'  m  ' 


(30) 


The  relative  uncertainty  in  the  imaginary  part  is  infinite  at  Af  =  0  because  the  mean  vanishes  at 
Af  =  0. 

D.  Evaluation  of  Variance  with  Noise  Present 

The  evaluation  of  the  variance  proceeds  along  the  same  lines  as  before  except  that  F^t)  has 
an  additional  noise  term  as  follows. 
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g  j  2 

F+(t)  =  H(f„  -  ^r,  t)  +  D(B  -  Af)  \  m  H(f  +  -  Af)  eZnivt  d v 

02  m  J-(B  /2)+Af  0  2 

m 


p-(Af/2)  +  (Bm/2) 
-(Af/2)-(B  /2) 


+  \  N(y  +  f  )  e2lril<t  dp 

%)_  (  a  /  o  \  I  rj  /  o  \  ” 


Af  r(B/2)-Af 

F"(t)  =  H(fn  +  -fi.t)  +  D(B  -Af)  \  m 

3  0  2  m  J-(Bm/2) 


tli/p  Af  .  a  Znivt  , 
H(f0 - y  +  Af)  e  dt> 


r(Af/2)+(B  /2)  ,  .  . 

+  \  m  N(,  +  f.)  e2™*  d,  . 

J(Af/2)-(Bm/2) 


(31) 


For  the  same  reasons  as  before  it  will  only  be  the  middle  term  on  the  right-hand  side  of  Eq.  (20) 
which  contributes.  In  the  expression  corresponding  to  Eq.  (22)  will  be  the  following  additional 
term  due  to  the  noise: 


sin7rB  (t  —  t’)  ,  * 

N  - r.~  -7j -  =  <F  (t)  F  (t’)> 

o  7r(t  —  t 1 )  n  n  a  vg 

Similarly,  there  will  be  an  identical  contribution  to  Eq.  (24),  namely. 


sin7rB  (t  —  t1)  * 

No  -  r(t-t')  -  <Fn'  (t,)  Fn(t»avg 


(32a) 


(32b) 


There  will  now  be  the  following  five  terms  in  expression  (25),  in  addition  to  those  already 
present: 


N 


sin7rBm(t  —  t’)l  2 
7r(t  -  t ' ) 


+  N 


sin7rB  (t  —  t') 
 m 


+  D(B  -Af) 


7r(t  -  t') 

MB  /Z)  -Af  9  ...  . 

(  nr  ,  —  27r  i(  t 1  -t)v 


R(0,  t  —  t’)  +  R(0,  t'  -  t) 


„  --  \  e 

m  J-(B  /2) 


<|H(f0,  v  +  Af)  |  2> 


avg 


nr 

Bm/2 


+  D(Bm  -  Af)  J  m' 


(B  /2)+Af 


dv  e2,rl(t'  <|  H(f0>  v  —  Af)  |  2> 


avg 


(33) 


With  the  assumption  [Eq.  (26)]  about  the  channel  properties,  these  terms  contribute  the  following 
amount  to  the  variance  after  the  double  integration  over  t  and  t’: 


B  N  0  +  2B  aN  0  +  D(B  -  Af)  2(B  -  Af)  a N  0 

mo  mo  m  m  o 

For  the  total  variance  in  the  presence  of  noise  one  therefore  obtains: 

2 


0B 


aZ  +  N2  +  2a N  +  D(B  -  Af) 
o  o  m 


(-r) 

'  m  / 


3a2  + 


2  a 


1  +  (27TT  Af) 


9  +  2a?N 
Z  o 


(34) 


The  signal-to-noise  ratio,  defined  as  the  square  of  the  mean  (excluding  the  noise  bias  which  is 
assumed  to  be  taken  out),  divided  by  the  variance  just  found  becomes: 
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(1)  For  the  real  part. 


(mean  signal)2  0  (^s\Z  /  _J_  \2 /sin7rTBm 

noise  (=  variance)  “  2Bm  \nJ  \1+{ZnT  ^)  \  *TBm  , 

l  m 

..LPs\2  1  i,laT.^  Fs  ll'1 

U2V  B2  i.Urr  A,)2  2NoBJ  ’ 


(35) 


where  ZaB  =  Pg  =  total  received  signal  power. 

(2)  For  the  imaginary  part,  the  corresponding  result  is  obtained  by 
multiplication  of  Eq.  (35)  by  (Zttt QAi)^ . 


When  the  signal  power  is  much  less  than  the  noise  power  within  the  bandwidth  Bm,  the 
result  [Eq.  (35)]  simplifies  quite  considerably  to  give 


(mean  signal) 
noise 


2 


0 

2B 


m 


_ 1 

1  +  (27rTQAf) 


sinffTB  \ 
- m  ] 

ttTB  / 


(36) 


It  is  interesting  to  observe  that  this  closely  resembles  the  signal-to-noise  ratio  of  the  output  of 
Price's  optimum  processor  in  the  case  where  the  signal  waveform  spans  many  coherence 

g 

intervals  (Price  ).  A  naive  interpretation  of  the  result  is 

2 

(signal)  _  J_  (signal  energy)  X  (signal-power  density) 

noise  2  ]  !  ~  \Z 

(noise-power  density) 


_1  g  0  /  signal-power  density ,  2 
2  ^m  '  noise-power  density 


(37) 


E.  Uncertainties  in  Target -Scattering  Function 


It  is  not  always  sufficient  to  be  able  to  discuss  the  variance  of  the  outcome  of  the  experimen¬ 
tal  scheme  shown  in  Fig.  1.  In  order  to  compare  matters  with  the  results  of  pulse  experiments 
designed  to  measure  the  target-scattering  function  <t(t,  f^)  directly,  it  is  necessary  to  find  the 
uncertainties  arising  in  the  determination  of  ct(t,  fp)  from  the  correlation  function  measurements. 
Suppose,  therefore,  that  we  have  gone  through  the  correlation  function  measurements  and  have 
determined  a  set  of  complex  numbers: 


V(mAf0,  ITq) 


m  =  0,  1,  .  .  . ,  M 
i  =  0,  1,  .  .  .  ,  L 


(38) 


and  that  the  integration  times  0  =  0^.  Our  estimate  of  the  target-scattering  function  obviously 
would  take  the  form 


<t(t,  f 


D^est 


2  0(m) 

f=-L 


M 

z 

m  =  -M 


27ri(mAfnT+iT  fn) 

V(mAf0,iTo)  e 


(39) 
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We  shall  not  be  concerned  with  scale  factors  here.  A  target-scattering  function  determined  in 
this  manner  will  deviate  from  the  true  scattering  function  for  two  reasons.  There  will  be  a 
truncation  error  because  the  summations  are  carried  to  finite  limits;  this  indicates  a  range 
resolution  of  the  order  At  =  f/MAfQ  and  a  frequency  resolution  Af^  =  f/LTQ  =  fjVL.  Also, 
there  will  be  a  systematic  error  due  to  the  finite  sampling  on  the  correlation  function.  Here 
we  shall  not  consider  these  systematic  errors.  Another  source  of  errors  is  caused  by  the  fact 
that  the  correlation  function  is  subjected  to  statistical  fluctuations  due  to  additive  noise  or  self¬ 
noise.  We  shall  be  concerned  with  these  latter  type  errors.  In  other  words,  we  would  like  to 
study  the  quantity 


M 


A(Ht,  fD) 


I  *  To}  -  <V(mAf0’ 1  To»avg> 


27ri(mAfr.r+i  T  f_.) 

U  O  JJ 


(40) 


-L  -M 


Hence  we  shall  attempt  to  determine  <Ac r(r,  f^)  ^avg-  It  is  found  that 


<A£r(r>fD)  >avg 


M 

l  ,  , 

m=-M  0(m)  I  n 


^  H  «V(mAf0,iTo)  V*(mAf0,  nTQ)> 
2iriT  f_(l-n) 


avg 


-  <V(mAf0'iTo)>avg  <V*(mAf0'nTo)>avg>  e 


o  D' 


(41) 


Assuming  here,  as  throughout,  that  the  channel  and  the  noise  are  both  Gaussian,  there  will  only 
be  the  terms  depending  on 

<  "O'  F3m|l'»  avg  <P3m«  *  T>  P;*  (f  t  T')>  avf  , 

which  will  contribute  to  the  entire  triple  sum  in  Eq.  (41).  The  first  of  the  two  factors  above  will 
be  identical  to  the  one  determined  previously  [see  Eqs.  (22)  and  (32a)]: 

27ritAf.m 


<F3m(t)  ^m^avg^’t'-^e 


D( B  -  mAfJ 
m  u 


>b  /2 

X  \  dy  e 


1 


27rii'(t’  -t) 


(B  /  2)+iAf. 
nV  0 


x<H*(fo-"T^'^-mAfo)H(fi 

0>f 


mAf^ 
0  2 


,  v  —  mAf 


o)>a 


-27rit'Af  m 

+  e  D(B  -  mAf. 

m  1 


B  /  2 
nv 


'Va,“o 


du  e 


avg 
27rii'(t'  -t) 


B  /2 


avg 


m 


+  D(Bm  -  mAfQ)  J  *“  dv  e 


27r  i^  ( t  *  -t) 


-(Bm/2)+fAf0 


sin7rB  (t  —  t1) 
m 


X  <H*(f0.  ,  -  mAf0)  H(f0,  „  -  mAf0)>avg  +  Nq  7r(t“\,) 


(42) 
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The  second  factor  takes  on  a  slightly  different  form  than  before.  Putting  T1  —  T  =  <5T  =  (f  —  n)  Tq, 
one  obtains 


<F3m(t)  F3m^'  "  6T»avg  =  R(°>  t  -  f  -  6T)  +  e 


-27ri(t,+6T)Afnm 


E»(Bm  -  mAf0) 


p(Bm/2)-mAf  ,  m^i0  v  /  im>i0  W 

J.(B  /2)  '  ^  A  0  +  ~~2  •  U  +  “*'<>)  Hif0  "  T~  ’  "  +  h^o)/ 


mAfn 


mAL 


avg 


X  e2^(t'-t+«T)  +  e27ritAf0m  D(B 

m  1 


0>  1 


(mAL  \  /  mAL  v 

f0  -  -T9  •  "  +  “*fo)  H(f0 +  -r2  •  *  +  mAfo)>avg 

r 


mAL 


dv 


^Tri^U’  +  ST-t) 


+  D(B  —  mAL)  \ 

m  °  J-<B  /2) 


(Bm/2|'mA'°  d.  tmAf0)|2>av 


+  N 


m' 


sin7rB  (t  —  t'  —  6T) 
_ m _ 

7r( t  -  t'  -  6 T) 


(43) 


In  the  subsequent  multiplication  and  integration  over  t  and  t',  only  those  terms  will  contribute 
which  depend  on  t  and  V  through  the  combination  t  —  t'.  This  leaves  the  following  eleven  terms 
(note,  we  put  t*  —  t  =  At): 


(I)  R(0,  At)  R*(0,  At  +  <5T)  + 

(n)  +  R(0,  At)  D(Bm  -  mAf0)  ^ 


p(B  /2)-mAf0  -2jrio(At+5T)  ,  , 

•  R(0,  v  +  mAL)  e  di<  + 


pB  /2 

(HI)  +  R*(0,  At  +  6T)  D(Bm  -  mAf0)  \  m 


-(B  /  2) 
nr 


(Bm/2)+mAf0 


R(0,  v  —  mAfg)  e27ri*'^'t  + 


-27rimAfn(At+6T)  pB  /2 

(IV)  +  D(B  —  mAL)  e  0  (  m 

m  o 


27riyAt 


\  R(mAfn,  v  —  mAL)  e  dv 

J-(Bm/2)+mAfo 


rl 

xi 


(Brn/2)-mAfn 


(V)  +  D(Bm  -  mAf0)  e 


’(B  /2)-mAfn 


m'  '  "  "0  Ar  ,  A  -  \  27rii^(At+6T)  , 

R*(mAf0,  v  +  mAfQ)  e  dy  + 

-27rimAfnAt  pB  /2 
u  i  m 


(B  /2) 

m 


1  m  R*(mAL,  v  —  mAf.)  e 

J-(B  /2)+mAfn 


27rii>At 


dy 


f|X3m/fc|"11I“A0  A.  .  Af  v  -27rii/(At+6T)  . 

\  R(mAfn,  v  +  mAL)  e  dv  + 

J-(B  /2) 


m 
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r*J  2 


(VI)  +  D(B  -  mAf  )  \  m' 

m  0  J-(Bm/2)+mAf0 


R(0,  v  —  mAfQ)  ^ 


X  ('<Bm/2)'mAf0  R(0j„  +  n^f  )  e-27rir(At+6T)  ^  + 


1 


(B  /2) 
m7 


sin7rB  At 

(vn)  +  n  - -Jp— 

o  7rAt 


R*(0,  At  +  <5T)  + 


sin7rB  (At  +  6T) 

(Vm)  +  No  — 6T)~  R(°’At)  + 


(EX)  +  D(Bm  -  mAf.) 

m  u 


N 


sin7rB  At  r*(B  /2)-mAfn 
m  (  m7  n 


- \  R(0,  i/  +  mAfn)  e 

°  *At  J-(B  /2)  0 

m7 

sin7rB  (At  +  <5T)  r*B  /Z 


-27riy(At+<5T) 


dv  + 


1  7T  +  0  1)  r 

™  *  "o  .(M  4  «T)  i 


m' 


'-(Bm/2)+mAf0 


R(0,  v  -  mAf0)  eZlTlvAt  di> 


9  sin7rB  At  •  sin7rB  (At  +  <5 T) 
<XI>  - "a,  ■  ,<At  4™T) - 


(44) 


Again  we  make  the  assumption  that  the  observation  time  spent  at  each  frequency  separation  is 
large  enough  so  that  the  integration  over  At  can  be  carried  between  infinite  limits.  Also,  in 
order  to  obtain  a  reasonably  simple  result  suitable  for  discussion  we  evaluate  all  the  above 
integrals  on  the  assumption  that  the  channel  properties  are  of  the  form  of  Eqs.  (26)  or  (27).  It 
is  found  that: 


Terms  (I)  +  (VII)  +  (VIII)  +  (XI) 

sin7rB  6T 


=  © 


™  —  (a2  +  2qN  +N2) 
7t6T  o  o 


Terms  (II)  +  (III)  +  (VI)  +  (IX)  +  (X) 

sin7r(B  —  mAf  )  6T 


=  ©D(B  -mAfJ 
m  0 


0' 


7t<5T 


(3a  +  2a N  )  ; 


Terms  (IV)  +  (V) 


=  2©D(B  -  mAf.) 

m  0 


sin7r(B  -  mAL)  <5T 
_ m _ u _ 

7t<5T 


1  +  (2t rr  Af) 
o 


(45) 


When  this  is  substituted  back  into  the  expression  for  the  variance,  we  obtain 
o  _  sin7rB  T  (f  —  n) 

<A^T.fD)  >avg~  Z  e(ST)  Z  Z 

-M  l  n 


,T  «  -  n)  <“  *  No>  *  DIBm  "  "“"V 


sin7r(B  -  mAL)  T  (i  -  n)  ? 

m  0  o  2  .  ^ 

X  - ^ -  (3a  +  2aN  + 


2or 


7rTQ(i  -  n) 


°  1  +  (27TToAf0f)2 


27riTQfD(f-n) 


(46) 
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Certain  conclusions  can  be  drawn  about  integration  times  0(m)  from  this  expression.  In  the 
general  case,  the  double  sum  over  I  and  n  is  a  function  of  both  m  and  f^;  let  this  function  be 
denoted  by  S(m,  f^).  Hence, 


<AcKt,  fD)  >avg 


M 

~  E 

-M 


S(m,  fD) 
0(m) 


(47) 


If  we  were  to  ask  which  set  of  0(m)  would  lead  to  an  optimization  of  the  output  signal -to-noise 
ratio  for  a  given  point  (t,  f^)  and  for  a  given  total  observation  time  0Q  =  2  0(m),  we  would 
obtain  m 


0(m)  -  Js( m,  fD) 


(48) 


However,  if  the  channel  is  distinctly  underspread,  or  if  the  noise -power  density  Nq  is  considerably 
greater  than  the  signal-power  density,  then  S(m,  f^)  will  be  essentially  independent  of  m  and  we 
conclude  that  equal  times  should  be  spent  at  each  frequency  separation. 


IV.  DISCUSSION  OF  APPLICATION  TO  RADAR  ASTRONOMY 


Let  us  finally  very  briefly  consider  an  example  of  the  determination  of  the  target-scattering 
function  in  radar  astronomy  by  this  method.  For  simplicity  we  consider  the  particular  case  of 
a  <  Nq,  i.e.,  where  the  energy  density  of  the  received  signal  is  less  than  the  noise-energy  density. 
This  is  a  situation  most  frequently  encountered  in  radar  astronomy  except  in  cases  of  lunar  and 
Venus  reflections  by  the  very  largest  radar  astronomy  installations.  It  is  also  a  situation  on 

g 

which  extensive  theoretical  considerations  regarding  optimum  detection  have  been  based.  In 

2 

the  variance  terms  in  Eq.  (46),  we  then  only  need  include  the  term  involving  N^  .  First  we  settle 
on  a  depth  range  to  be  investigated;  ideally,  this  depth  should  be  equal  to  the  radius  of  the  target. 
In  practice,  however,  we  might  specify  the  depth  range  somewhat  smaller  than  this  because  the 

echo  toward  the  limbs  of  the  object  is  frequently  too  weak  to  be  detected  anyway.  We  let  the 

-1 

depth  range  be  Afn  ;  this  specifies  the  smallest  frequency  separation  we  want  to  use  in  the 
experiment.  Similarly,  the  smallest  time  shift  Tq  must  be  chosen  about  equal  to  Bm  .  The 
next  question  is  concerned  with  the  number  of  resolution  cells  we  want  within  the  specified  am¬ 
biguity  range.  Suppose  we  desire  2M  resolution  cells  in  depth  and  2L  in  frequency.  With  the 
flat  noise  spectrum  assumed  in  the  above  calculations  we  obtain 


,  A  /  ,  \2\  m2  2M  +  1 

<  Acr(r,  fD)  >ayg  =  No 


0 


L  L 

E  E 

f=- L  n=-L 


sinird  —  n) 

T  v(t  -  n) 


2?riT  f_(l-n) 
o  u 


(49) 


_  \ 

If  L  is  at  least  moderately  large,  the  double  sum  over  l  and  n  will  be  approximately  (2L  +  1)  Tq 
when  |fp|  <  Bm/2,  and  zero  otherwise.  Hence,  within  the  bandwidth  containing  the  signal  we 
have 


<Acr(-r,fD)2>avg  =  N2(ZL  +  1)  (2M  +  L) 


(50) 


On  the  assumption  of  an  exponential  target-scattering  function,  the  mean  value  of  the  signal  will 
be 


<  cr>  avg  -  #  AE7 


-t/t 


B 


when  f 


0  o 


D1 


m 

2 


(51) 
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Now,  Af^r^  is  some  numerical  factor  depending  on  the  ratio  of  the  decay  constant  r  and  the 
depth  to  which  we  want  to  investigate  the  target.  Let  this  ratio  be  (3  The  signal-to-noise 
ratio  becomes 


signal  2  = 
'  noise  ' 


Bm0^ 


(2L  +  1)  (2M  +  1) 


e 


-2  t/t 


o 


(52) 


where  (2L  +  1)  (2M  +  1)  is  the  number  of  resolution  cells.  Finally,  introducing  the  total  observa¬ 
tion  time  =  M0  and  the  total  power  P  =  aB  2  we  obtain 
o  ^  rec  m 


signal  2  _  /_L  \2 
'  noise  '  \2No/ 


/J2e 


-2t/t 


M(2L  +  1)  (2M  +  1)  B 


(53) 


V.  CONCLUSIONS 

We  have  shown  that  a  measuring  scheme  can  be  devised  which  makes  the  two-frequency 
method  of  investigating  statistically  stationary  spread  channels  possible  both  in  the  underspread 
and  overspread  situations.  In  case  the  channel  has  correlated  taps,  the  statistical  properties  of 
the  channel  vary  with  center  frequency  fQ  so  that  the  observation  must  be  repeated  for  several 
center  frequencies.  The  method  described  appears  to  be  a  rather  flexible  one  in  that  the  various 
channels  in  the  observational  scheme  shown  in  Fig.  1  should  be  rather  easy  to  implement  and  to 
alter  at  will,  particularly  if  some  of  the  processing  is  being  carried  out  by  digital  techniques. 

It  appears  to  be  rather  difficult  to  compare  the  method  with  other  methods  of  obtaining  target¬ 
scattering  function  estimates  because  few  such  estimates  are  available  and  the  constraints  on 
the  measuring  methods  are  often  different.  One  might  think  that  a  direct  comparison  could  be 

4 

made  of  the  results  of  Gallager's  chirp  experiment  and  the  present  results.  It  turns  out, 
however,  that  it  is  not  at  all  obvious  how  such  a  comparison  is  to  be  carried  out  because  in  the 
present  report  we  measure  R(Af,  T)  for  a  particular  Af  with  variable  T,  whereas  the  chirp 
method  measures  R(aT,  T)  where  a  is  a  constant  determining  the  rate  of  change  of  the  fre¬ 
quency  in  the  chirp.  We  believe,  however,  that  the  method  described  might  be  advantageous 
in  cases  where  there  is  a  peak-power  constraint  on  the  transmissions  that  can  be  used.  A  dis¬ 
advantage  under  these  circumstances  would  be  the  fact  that  the  average  transmitted  power  is 
only  equal  to  half  the  peak  power. 
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